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We study the transport properties of a hnite three dimen- 
sional disordered conductor, for both weak and strong scat- 
tering on impurities, employing the real-space Green func- 
tion technique and related Landauer-type formula. The dirty 
metal is described by a nearest neighbor tight-binding Hamil- 
tonian with a single s-orbital per site and random on-site 
potential (Anderson model). We compute exactly the zero- 
temperature conductance of a finite-size sample placed be- 
tween two semi-infinite disorder-free leads. The resistivity is 
found from the coefficient of linear scaling of the disorder- 
averaged resistance with sample length. This "quantum" re- 
sistivity is compared to the semiclassical Boltzmann expres- 
sion computed in both Born approximation and multiple scat- 
tering approximation. 

PACS numbers: 72.15.Eb, 72.15.Lh, 72.15.Rn 



Ever since Anderson's seminal paper ,0 a prime model 
for the theories of the disorder induced metal-insulator, 
or localization-delocalizationB (LD), transition in non- 
interacting electron systems has been the tight-binding 
Hamiltonian (TBH) on a hypercubic lattice 
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with nearest neighbor hopping matrix element t between 
s-orbitals (r|m) = ip(r — m) on adjacent atoms located at 
sites m of the lattice. The disorder is simulated by tak- 
ing random on-site potential such that e m is uniformly 
distributed in the interval [-W/2,W/2]. This is com- 
monly called the..." Anderson model" . There are many 
numerical studiesu of the LD transition, which occurs in 
three-dimensions (3EL) for a half-filled band at the criti- 
cal disorder strengths W c w 16. 5t. Experiments on real 
metals with strong scattering or strong correlations of- 
ten yield resistivities which are hard to analyze. Theory 
gives guidance in two extreme regimes: (a) the semiclas- 
sical case where quasiparticles with definite k vector jus- 
tify a Boltzmann approach and "weak localization" (WL) 
correction!] and (b) a scaling regimeO near the LD transi- 
tion to "strong localization" . Lacking a complete theory 
it is often assumed that the two limits join smoothly with 
nothing between. Experiments, however, are very often 
in neither extreme limit. The middle is wide and needs 
more attention. 

Here we give a 3D numerical analysis focused not on 
the transition itself but instead on the resistivity for 
1 < W/t < W c /t; specifically we ask how rapidly does the 
resistivity p(W) deviate from the values predicted by the 
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FIG. 1. Resistivity p at Ef = (lower panel) and 
Ep = 2 At (middle panel), from a sample of cross section 
A — 225 a 2 , normalized to the semiclassical Boltzmann resis- 
tivity pb calculated in the Born approximation. Also plotted 
are the ratios of pe to the Boltzmann resistivity px obtained 
using a T-matrix for multiple scattering on a single impurity. 
The upper panel shows putative mean free paths obtained 
from pa (labeled by B) or px (labeled by T). Error bars at 
small W/t are smaller than the size of the dot. 



usual Boltzmann theory valid when W It has long 
been assumed that "Ioffe-Regel condition" u£ ~ 1/kp ~ a 
(£ being the mean free path, and a being the lattice con- 
stant) gives the criterion for sufficient disorder to drive 
the metal into an Anderson insulator. Figure ^ shows 
that this is wrong. By W/t ~ 4, where I is close to a, 
there is little sign of a divergence away from the semiclas- 
sical extrapolation, and the LD transition is postponed 
to much larger values of W/t. A cleaner discussion is pos- 
sible using Kubo theory, which does not define £, but al- 
lows a definition of the diffusivity Di of an eigenstate \i), 
as shown below in Eq. (|J). In the semiclassical regime, 
Di — > Dk = The diffusivity Dk diminishes as 

(W/t)~ 2 in Boltzmann theory. As l/a approaches a min- 
imum value (~ 1), Di decreases toward D m - ln = ta 2 /h, 
which can be regarded as a minimum metallic diffusiv- 
ity below which localization sets in. But there is a wide 
range of W/t over which Di < D min and yet the Boltz- 
mann scaling D ~ (t/W) 2 is approximately right. In this 
regime single particle eigenstates |i) are neither ballisti- 



1 



cally propagating nor are they localized. There is a third 
category: intrinsically diffusive £3 A wave packet built 
from such states has zero range of ballistic motion but 
an infinite range of diffusive propagation. Such states 
are found not only in a narrow crossover regime but over 
a wide range of parameters physically accessible in real 
materials and mathematically accessible in models like 
the Anderson model. In this regime, there is not a sim- 
ple scaling parameter nor a universal behavior. But the 
behavior is quite insensitive to a changes in Fermi energy 
Ep or kgT, and scales smoothly with W/t. 

The traditional tool for computation of p has been the 
Kubo formula,H originally derived for a system in the 
thermodynamic limit. In a basis of exact single particle 
electron state \i) of energy a, this can be written as 

where is the sample volume, / is the equilibrium Fermi- 
Dirac distribution, N(Ef) the density of states at Ef, 
D the mean diffusivity, and state diffusivity is given by 

Di = irh^2\(i\v x \j)\ 2 5(ci - ej) (3) 

3 

where v is the velocity operator. These formulas, while 
correct, are hard to use numericallyB_Thanks to the re- 
cent advances in mesoscopic physics, E3 it is, now appar- 
ent that the Landauer scattering approacblij (or equiva- 
lent "mespscopic" Kubo reformulation for the finite-size 
systemsEij) provides superior numerical efficiency when 
computing the transport properties of finite disordered 
conductors. It relates the conductance of a sample to 
its quantum- mechanical transmission properties. This 
formalism emphasizes the importance of taking into ac- 
count the ipterfaces between the sample and the rest of 
the circuit .Lj Transport in the sample is phase-coherent 
(i.e. effectively occurring at zero temperature); the dis- 
sipation and thus thermalization of electrons (necessary 
for the establishment of steady state) takes place in other 
parts of the circuit. 

Our principal result for the (quantum) resistivity 
of Anderson model, using Landauer-type approach, is 
shown on Fig. |] for two different Fermi energies Ep = 
(half-filled band) and Ep = 2At (approximately 70% 
filled band but falling somewhat as W, and thus the 
band- width, increases). The linearized Boltzmann equa- 
tion — eE-Vk<9//<9ek = (dFk/dt) S catt serves as a reference 
theory. Here ek is the energy band for W — 0, namely 
e k = 2t^cosfc Q , hvk a is de/dk a , and Fk is the non- 
equilibrium distribution. The collision integral is 




= -irEl^'| 2 ( Fk -^'Wek-e k ,). (4) 



The mean squared matrix element of the random po- 
tential |Vkk'| 2 i m Born approximation, is = W 2 /12, 



where (...) denotes average over probability distribu- 
tion P(e m ) = (l/W)6(W/2 - |e m |). This equation as- 
sumes that quasiparticles propagate with mean free path 
I > a between isolated collision events. The equation 
is exactly solvable, yielding (for fc^T <C t) 1/pb = 
e 2 T(n/m) cS , with (n/m) cS = J2 v lx^( e k - E F )/£l, and 
H/t = 2ttN(E f )W 2 i '12. We have evaluated (n/m) cS 
and N(Ep) numerically. To within factors of order one, 
the Boltzmann-Born answer for the semiclassical resis- 
tivity is p B = (nha/e 2 ){W/U) 2 . When W = 3t and 
a = 3A, Pb is 125 pQcm, typical of dirty transition metal 
alloys, and close to the largest resistivity normally seen in 
dirty "good" metals. Figure |l| plots p/pb versus (W/t) 2 . 
Even for W = lOt there is less than a factor of 2 devia- 
tion from the (unwarranted) extrapolation of the Boltz- 
mann theory into the regime W > t. Boltzmann theory 
can be "improved" by including multiple scattering from 
single impurities, that is, replacing the impurity poten- 
tial by the T-matrix T m (z) = e m /(l — s m g(z)) where 
g(z) — (l/N s )Y2(z — £k) _1 is the free particle Green 
function (N s is the number of lattice sites) . To next or- 
der the mean square T-matrix is 

W 2 ( 3W 2 \ 

|T m (z)P = V — ( 1 + —(gg* + gg + g*g*) + . . . J , 

(5) 

where the first term is the Born approximation and the 
coefficient of the correction (~ 0(W A )) changes sign 
from negative to positive as Ep moves from to 2 At. 
As shown on Fig. the resistivity does not behave like 
|T m (,z)| 2 ; multiple scattering with interference from pairs 
of impurities is at least equally important, and the "ex- 
act" p(W) is less sensitive to details like Ep than is the 
T-matrix approximation. The rest of the paper presents 
the method of calculation and describes a bit of meso- 
scopic physics of very dirty metals. 

The central linear transport quantity in the meso- 
scopic viewjj as well as in the scaling theory of lo- 
calization^ is conductance G rather than conductivity 
<j(L) = L 2 ~ d G(L) (the bulk conductivity is an inten- 
sive material constant defined only in the thermodynamic 
limit, a = limL->oo L 2 ~ d G(L)). We use a Landauer-type 
formula to get the exact quantum conductance G of fi- 
nite samples with disorder configurations chosen by a 
random number generator. Finite-size samples permit 
exact solutions for iaip-| strength of disorder. Similar to 
other recent works,ESE3 the bulk resistivity is extracted 
from the disorder-averaged resistance (R) by finding the 
linear (Ohmic) scaling of (R) versus the length of the 
sample-.L at fixed cross section A (Fig. ||). Two kinds of 
errorsEj may arise: (a) The transition from the Ohmic 
regime to the localized regime occurs for length of the 
sample L ~ £ which happens whentll G = 0(2e 2 /h). If 
L is made large enough, G will always diminish to this 
magnitude. Therefore, we avoid using the sample sizes 
with too small G. (b) Fi»ite-size boundary conditions 
and non-specular reflectionta cause density of states and 
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scattering properties of the sample to be slightly altered 
as compared to the true bulk. We expect these effects 
to be small for our samples where £ is smaller than the 
transverse size vA 

A two probe measuring configuration is used for com- 
putation. The sample is placed between two disorder-free 
(e m = 0) semi-infinite leads connected to macroscopic 
reservoirs which inject thermalized electrons at electro- 
chemical potential p,L (from the left) or (from the 
right) into the system. The electrochemical potential dif- 
ference eV = (J,l— hr is measured between the reservoirs. 
The leads have the same cross section as the sample. The 
hopping parameter in the lead and the one which couples 
the lead to the sample are equal to the hopping parame- 
ter in the sample. Thus, extra scattering (and resistance) 
at the sample-lead interface is avoided but transport at 
Fermi energies \Ep\ greater than.-the clean-metal band 
edge | | = 6t cannot be studied El Hard wall boundary 
conditions are used in the y and z directions. The sam- 
ple is modeled on a cubic lattice with N x N y x N z sites, 
where N y = N z = 15 and lengths L = Na are taken from 
the set N e {5,10,15,20}. 

The linear conductance is calculated using an expres- 
sion obtained from the Keldysh technique^ 



4e 2 
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Here ImEx^, = (J7 L R — #)/2i are self-energy matri- 
ces (r-retarded, a-advanced) which describe the coupling 
of the sample to the leads, and G\ N , G a N1 are Green 
function matrices connecting the layer 1 and N of the 
sample: G r ' a = (E — H — E^ )" 1 (G Q = [G r ] t ), with 
£ r = ± r L + ± r R (± a = [S r ]t). The self-energy matrices 
introduced by the leads are non-zero only on the end laj^-. 
ers of the sample adjacent to the leads. They are givenE2l 
by S£ )fl (n, m) = t 2 g r LR (n s , ma) with g r LtR (n s , m s ) be- 
ing the surface Green functionEj of the bare semi-infinite 
lead between the sites ns and mj in the end atomic layer 
of the lead (adjacent to the corresponding sites n and m 
inside the conductor). Positive definiteness of the oper- 
ators — 2ImSi /j R makes it possible to find their square 
root and recast the expression under the trace of Eq. (||) 
as a Hermitian operator. The expression (^) then looks 
like the Landauer formula involving the transmission ma- 
trix t 
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t = 2V-ImE L G^V" ImE 
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(8) 



or transmission eigenvalues T n when the trace is evalu- 
ated in a basis which diagonalizes tt^. 

For the case of two probe geometry the average trans- 
mission in the semiciassical transport regime (a < £ <C 
L < f) is given byM (T) = £ /(£ Q + L), with £ being 
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FIG. 2. Linear fit (R) = Ci + pL/A, (A = 225 a 2 ) for the 
disorder-averaged resistance (R) in the band center Ef = 
and different disorder strengths W . The intercept Ci is de- 
creasing with increasing W (i.e. it is not determined just by 
the contact resistance Tvh/147e 2 ) and becomes negative for 
W > 7t. The inset shows examples of the distribution of re- 
sistances Pl(R) (for L = 15a) versus logi?. The distribution 
broadens either by increasing W or the length of the sample 
(the units on the y-axis are arbitrary and different for each 
distribution) . 



of the order of I. Thus, the semiclassical limited of the 
Landauer formula for conductance (G) — (e 2 /irh) M (T) 
(measured between points deep inside the reservoirs) in 
the case of not too strong scattering should have the form 



(GY 



(9) 



It describes the (classical Vseries addition of two resistors. 
The "contact" resistancetil Rq = nh/e 2 M is non-zero, 
even in the case of ballistic transport when the second 
term containing the resistivity p = (Trh/e 2 ) A/£qM van- 
ishes. Here M ~ kpA is the number of propagating 
transverse modes at Ep, also referred to as "channels". 
A ballistic conductor with a finite cross section can carry 
only finite currents (the voltage drop occurs at the lead- 
reservoir interface). Using this simple analysis for guid- 
ance, we plot average resistances (taken over iV con f = 200 
realization of disorder) versus L in Fig. g, and fit with 
the linear function 



(R) = C x + C 2 L, C 2 = p/A. 



(10) 



The resistivity p on Fig. [I] is obtained from the fitted 
value of C 2 ■ For very small values of W the constant C\ is 
approximately equal to Rc — Trh/e 1 M (where M = 147 
is the number of open channels in the band center) . To 
our surprise, C\ diminishes steadily with increasing W, 
and even turns negative around W ^ It. 

The quantum conductance G fluctuates from sam- 
ple to-sample exhibiting universal conductance fluctu- 
ationsE 2 ] (UCF) AG = VVaTG ~ e 2 /irh in the semi- 
classical transport regime G ^ e 2 /nh. The inset on 
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FIG. 3. The conductance fluctuations (AG = VvarG at 
Ef = 0) from weak to strong scattering regime in the disor- 
dered cubic samples 10 x 10 x 10 and 15 x 15 x 15. 



It is interesting to note that in many d-band inter*, 
metallic compounds, p "saturates" at a constant valuet-j 
rather than following the semiclassical extrapolation, 
that is, increasing linearly with T at high T. High T c 
materials ancLdoped Ceo metals, on the other hand, do 
not saturateo Within Boltzmann theory, the static dis- 
order measured by (W/t) 2 plays the same role as thermal 
disorder or squared lattice displacement oc ksT. Our nu- 
merical results thus can be described as "failing to sat- 
urate." Similar failure was seen-in high T Monte Carlo 
studies by Gunnarsson and HanBa 

We thank I. L. Alciner for interesting discussions and 
challenging questions. Suggestions provided by J. A. 
Verges are acknowledged. This work was supported in 
part by NSF grant no. DMR 9725037. 



Fig. H shows the distribution of resistance^! Pl{R) for 
our numerically generated impurity ensemble. The error 
bars, used as weights in the fit (Eol), arc computed as 
S(R) = ^Va,rR/N con{ . We find that AG is indeed inde- 
pendent of the size L (of cubic samples), but decreases 
systematically by a factor w 3 as W increases to the crit- 
ical value W c (Fig. ||). On the other hand, AR, being 
similar to AG/G 2 , depends on sample size. As W ap- 
proaches W c , G gets smaller until (for our finite samples) 
AG/G ~ 1. At this point the distribution of resistances 
R = 1/G becomes very broad and (R) begins to rise 
above 1 /(G). For L = 15 this happens when W > 12t. 
At large W the conductance of long samples (N = 20) 
becomes close to e 2 /irh and deviations from Ohmic scal- 
ing are expected. Therefore, we do not use these points 
in the fitting procedure when W > lOi (keeping the con- 
ductanceEj of the fitted samples G > 2e 2 /irh). 

We do not have a complete explanation for the de- 
viation of C\ (10) from the quantum point contact re- 



sistance Rc- In the semiclassical regime G 3> e /irh 
there are corrections to the Ohmic scaling G oc L d ~ 2 . 
The Diffuson-Cooperon diagrammatic perturbation the- 
ory gives a (negative) WL correction^ 



a(L) = a 



1 



ir 2 h2V2L ir 3 h£' 



(11) 



where £' is a length of order I (its precise value does 
not lead to observable consequences in the experiments 
studying WL, as long as it is unaffected by the temper- 
ature and the magnetic field). The positive l/L term in 
Eq. [ll] provides a possible picture for our finding that 
Giin Eq. ( ^fj| ) goes negative as W increases. However, 
this picture is an extrapolation from the semiclassical 
into the "middle" regime of intrinsically diffusive states, 
and therefore should be given little weight. The negative 
values of G\ is better regarded as a new numerical result 
from the mesoscopic dirty metal theory. 
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